IAO-A08«  071 
UNCLASSIFIED 


PRINCETON  UN2V  NJ  OCPT  Of  STATISTICS  p/s  ,»/, 

A  ROBUST  CONFIDENCE  INTERVAL  FOR  SAMPLES  OF  FIVE  OBSERVATIONS. (U) 

*>v  ”  *  KiFAD**  DAAA2R-76-A-0298 

TR-15A-SCR-2  ARO-1A2M.1A-M  NL 


ARO-I*2**.l«-N 


AO  A088071 


UNILASS I F I L  D 


X.  ^  f  <  f  J 


SECURITY  CLASSIFICATION  OF  THIS  PACT  (*T>«n  Omtm  F.ntm rmd) 


a  r~n 


rf) 


tM, 


REPORT  DOCUMENTATION  PAGE 


'  READ  INS' 
fcBKORE  COI 


ORM 


'ORT  NUMBER 


IA2A4.  1 4-M 


1.  GOVT  ACCESSION  NCJ; 


3.  RfClPlEN'T’S  C  A^i 


TIE  (and  Subtltlm) 


5.  REPORT  4  PE 


EREO 


4  ROBUST  CONFIDENCE  INTERVAL  FOR  SAMPLES  OF 
J I VE  OBSERVATIONS  -  *  **  | 


Technical 


6.  PEfiE-ORMING  ORG.  REPORT, 


i2s 


■.;«qm4nact  or  grant  num 


u 


DAAG29-76-G-^98^  jN 


».  performing  organization  name  and  address 

Princeton  University 
Princeton,  NJ  085^0 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  A  WORK  UNIT  NUMBERS 


It.  CONTROLLING  OFFICE  NAME  AND  AOORESS 

U.  S.  Army  Research  Office 
Post  Office  Box  12211 
Research  Triangle  Park,  NC  27709 


.12*.  REPORT  DATE 


Nov  79 


U.  NUMBER  DT  PAGES 


ts/p 


14.  MONITORIN  J  JHT6NCY  NAME  •  ADDRESS (II  dllloront  from  Controlling  Oltlco) 

(/^avi 


IS.  SECURITY  CLASS,  (of  thlm  raporf) 

Unclassified 


4 — - 


I  5  a.  DECLASSI  Ft  CATION/ DOWN  GRADING 
SCHEDULE 


««.  DISTRIBUTION  STATEMENT  fof  t hit  Roport) 

Approved  for  public  release;  distribution  unlimited. 


.  ,i  e  W 


V  AUG  1  8  1980 


rp 


W 


17.  DISTRIBUTION  STATEMENT  (ot  tha  abstract  antarad  In  Block  20.  it  dltlarant  from  Rapori) 


NA 


IB.  SUPPLEMENTARY  NOTES 

The  view,  opinions,  and/or  findings  contained  in  this  report  are  those  of  the 
author (s)  and  should  not  be  construed  as  an  official  Department  of  the  Army 
position,  policy,  or  decision,  unless  so  designated  by  other  documentation. 


19.  KEY  WORDS  (Conilnua  on  old*  II  nacaaaary  and  idantity  oy  otocm  numttar) 


confidence 
sampl i ng 
stat i st i cs 
probabi 1 i ty 


simulation 

Monte  Carlo  methods 


ABSTRACT  fCmmtBmaa  am  mamma  «Mi  ft  wwify  amd  Idantity  by  bloc  k  mtmkar) 

A  robust  confidence  interval  using  biweights  for  the  case  of  five  observations  is 
proposed  when  the  underlying  distribution  has  somewhat  heavier  tails  than  the 
Gaussian.  The  distribution  of  a^^like  statistic  is  approximated  by  a  Student’s 
t  on  the  nominal  four  degrees  of  freedom  using  different  scale  factors  which 
depend  upon  the  value  of  the  blweight  weights.  Results  given  by  Monte  Carlo 
simulations  indicate  that,  even  for  very  high  coverage  probabilities,  the 

intervals  proposed  are  highly  efficient,  in  terms  of  the  expected  length  of  the 
confldencp  Interval. 

DO  •  JAM*  71  1473  cornokor  •  nov  » is  obsolete 


UNCLASSIFIED 

crniwrv  n  •ifiririTfna  or  tm 


80 


0 


i 


"r 


A  ROBUST  CONFIDENCE  INTERVAL 
FOR  SAMPLES  OF  FIVE  OBSERVATIONS 


by 


Karen  Kafadac 
Princeton  University 
and 

Oregon  State  University 


Technical  Report  No.  154,  Series  2 
Department  of  Statistics/ 
Princeton  University 
November  1979 


This  research  was  supported  in  part  by  a 
contract  with  the  U.  S.  Army  Research  Office, 
No.  DAA329-76-0293 ,  awarded  to  the  Department 
of  Statistics,  Princeton  University, 
Princeton,  New  Jersey. 


( 


ABSTRACT 


A  robust  confidence  interval  using  bi¬ 
weights  for  the  case  of  five  observations  is 
proposed  when  the  underlying  distribution  has 
somewhat  heavier  tails  than  the  Gaussian.  The 
distribution  of  a  "t"-like  statistic  is  approx- 

t 

imated  by  a  Student's  t  on  the  nominal  four 
degrees  of  freedom  using  different  scale  fac¬ 
tors  which  depend  upon  the  value  of  the  bi¬ 
weight  weights.  Results  given  by  Monte  Carlo 
simulations  indicate  that,  even  for  very  high 
coverage  probabilities,  the  intervals  proposed 
are  highly  efficient,  in  terms  of  the  expected 
length  of  the  confidence  interval. 
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Q.  INTRODUCTION 

In  an  earlier  report,  the  author  [9]  considered  the  use  of 
the  biweight  in  constructing  a  confidence  interval  for  a  sample 
of  at  least  ten  observations.  Using  the  Student's  t  critical 
point  on  nine-tenths  of  the  nominal  degrees  of  freedom,  it  was 
found  that  the  efficiency  of  a  100«(l-o)%  confidence  interval 
in  the  Gaussian  and  in  symmetric  stretched-tailed  situations  ex¬ 
ceeded  80%  across  a  wide  range  for  a.  In  this  report,  we  brave¬ 
ly  explore  the  performance  of  the  biweight  in  a  "t"-like  statis¬ 
tic  when  we  have  only  five  observations.  We  are  looking  for 
good  performance,  not  only  in  the  (unlikely)  event  that  our 
sample  is  truly  Gaussian,  but  also  if  our  sample  comes  from 
a  population  with  somewhat  heavier  tails  than  the  Gaussian. 

Little  is  known  about  the  results  of  robust  procedures 
of  the  location  problem  alone  on  such  small  6ize  samples. 

The  Princeton  Robustness  Study  [l3  concluded  that,  in  terms  of 
95%  confidence  intervals,  the  estimates  could  show  considerable 
differences  in  non-Gaussian  situations  (Section  7B);  their 
recommendation  was  a  redescending  Hampel-type  estimator  (Sec¬ 
tion  6L).  Huch  of  the  literature  on  the  interval  problem  for 
small  samples  has  concentrated  on  the  analytic  distribution  of 
Student's  t  statistic  (e.g.,  [1),  (73).  Tor  more  general 
stretched-tailed  situations,  several  authors  have  shown  that 
Student's  t  is  highly  conservative  (e.g.,  [13],  CIS]).  Except 
for  specific  underlying  densities,  a  general  solution  to  the 
interval  problem  has  not  been  considered.  The  situation  is 


particularly  complicated  by  the  fact6  that 

(i)  suspect  ‘'outliers,"  even  in  Gaussian  samples, 
are  not  uncommon; 

(ii)  a  95%  confidence  interval  for  five  Gaussian 
observations  necessarily  extends  beyond  the 
range  of  the  data; 

(iii)  an  extremely  heavy-tailed  situation  offers  Just 

minimal  amount  of  information  required  for  a  con¬ 
fidence  interval,  for,  although  the  variance  of  a 
given  H-estimate  is  finite,  higher  moments  may  not  be. 

While  there  exist  many  estimates  to  use  in  constructing  robust 
confidence  intervals,  this  report  considers  only  the  biweight  in  a 
"t"-like  statistic,  largely  on  the  basis  of  its  previous  success  in 
problems  of  interval  estimation  ([6],  (93),  regression  ((23),  and 
time  series  ((33).  This  report  is  divided  into  three  parts; 

Part  A  presents  the  results  of  biweight-"t"  in  the  three  sampling 
situations;  Part  B  investigates  a  method  to  improve  our  estimate 
of  the  variance  of  the  biweight  via  "compartmentalizing,"  and 
Part  C  offers  conclusions  and  strategies  for  the  case  of  five 
observations . 


ON  "UNUSUAL**  SAMPLES. 


PAKT  A. 


THE  FAILURE  OP  “t“ 


1 .  For m  of  biwe  ight-^t*  and  concepts. 


Foe  a  definition  of  the  biweight  and  its  associated 
variance,  the  reader  is  referred  to  ll2lj  we  mention  here 
only  the  computational  methods.  The  biweight  estimate  of 
location,  ^bi'  is  defined  as  the  solution  to  the  equation 
n 

3  +  ((*,  -  Tb.)/(CS)!  -  0  ,  (Li 

l-l 


where 


uU-uV 

0 


u*w(u)  , 


I  U  I  <  1 
else 


Here,  s  is  an  estimate  of  scale  from  the  sample  Xj,,***,  xR, 
and  c  la  a  multiple  of  the  scale.  ( K  choice  of  c  recom¬ 
mended  in  I12I  is  that  for  the  denominator,  c*s  ,  is 
between  4o  and  6a  in  the  Gaussian  case.  In  this  study 
we  will  choose  c  such  that  c*s  is  roughly  6a  for  the 
Gauss  ian. ) 

We  may  rewrite  (1)  in  terms  of  the  "weight  function”, 
w(u) ,  where 


w(u)  *  *(u)/u  , 


n 


bl 


S  x.wtu.l 
i-i 


(2> 


r  w(u.) 
i*l  1 

Equation  (2)  suggests  an  iterative  solution.  We  start  the 
iteration  with  a  robust  estimate  of  location  (in  this  study, 
the  median  of  the  sample).  The  location  estimate  at  the  kth 
iteration,  ,  k  >  1 ,  is  found  by 


5  «,«(<*,  - 

t(M  .  _ 

bl  n 

3  w((xt  -  T^"M)/(c*8|) 


Ol 


In  determining  an  estimate  of  scale  to  use  in  (3),  former 
studies  (see,  e.g . , [l  ]  ,[ll )  )  suggest  the  median  absolute 
deviation  from  the  median  (MAD)  » 


.(0) 


med  I  x . 


T(0) 

Tbi 


Por  reasons  to  become  clear  later.  Lax  (111  showed  that  a 
more  efficient  scale  estimate  may  be  that  using  the  func¬ 
tional  form 

=bi  -  «U2,«vl0’ri,m,(ai"  <«> 

where 


a 


l 


whence 


q4ll<<ai» 


5  +  (Oj) 
l-L 


[LHIr  i-i*j  *•<’.>>] 


Here,  as  before,  T*  is  the  median  of  the  sample,  s'^ 
is  the  M*D,  and  cQ  is  again  chosen  in  order  that  c0's<0> 
is  approx imately  the  desired  multiple  of  a  in  the  Gaussian 
case.  (Since  s*°*  Z  (2/3) a  for  a  Gaussian  sample,  we 
choose  cQ  -  9  for  this  calculation.) 

Finally,  the  denominator  of  our  "t  "bl  statistic  is 
given  by  ,  where  Sbl  estimates  the  variance  of  . 

Huber  (81  derives  the  theoretical  asymptotic  variance  of 
Tbl'  tcom  «hich  we  may  obtain  a  finite-sample  approximation 


sbl  "  vic<Tbi>  “  <c*«bir’(tai^» 


u4  -  , 

1  c  »bi 

M  in  e^uat Ion  M).  Notice  that,  in  funct  tonal  foci. 


5  (X  - 


.  claeaicai  aanpie  a* 


in  the  5ausa ian  case.  However,  s£(  uses  the  swdian  and  the 


location  and  scale  estiaates  Tbj  and  sbl>  Notice  also  that 
93  defined  in  (Si  may  be  written 


i  ufil-ufl9 


,iii,lui11  -  7^ 


••**■*.—  n  ~|(“  n  — « 

5  (1-Uj) (l-5u|)  eas(  1,-1+  5  (l-u|l Cl-Su |> 1  I 

The  exponents  of  (1-u*)  ,  (l-u*>  ,  and  (l-SoJ)  ,  respec¬ 

tively.  suggest  the  subscript  and  the  name  Mil  wtdther-  for 
*bi  (Equation  (|))4  our  blwe  ight-"t”  statistic  then  takes 


2.  Results  on  samples  of  size  five. 

Performance  of  biweight-M"  will  be  evaluated  on  three 
different  distributions: 

°  Gaussian 

°  One-Wild  (H  observations  from  N(0,l)i 

1  unidentified  observation  from  N(0.  100)) 

°  Slash  (N(0,1)  deviate  /  independent  U[0,1]  deviate)  . 
These  three  situations  are  likely  to  cover  a  reasonably  broad 
range  of  stretched-tailed  behavior.  The  critical  points  of  the 
distribution  were  all  computed  via  a  Monte  Carlo  swindle,  the 
results  of  which  may  be  found  in  [5].  There  were  6M0  samples 
in  the  simulation  for  each  sampling  situation. 


Var(X) 


The  success  of  biweight-Mt"  will  be  measured  primarily 


in  terms  of  "efficiency" 
length  (ECI1.),  i.e.  , 


et'f(a)  = 


of  the  expected  confidence  interval 


ECIL  .  (a) 
min 

ECIL  7  r?aT 
_  actual 


where  ECIL(u)  was  defined  by  Gross  ([2])  as 


ECiL(a)  =  2‘ai-point ’ave(denominator  of  "t") 


and  ECIL  .(a)  is  the  "shortest"  obtainable  for  the  situation 
min 

(see  (9)).  Furthermore,  we  shall  be  interested  in  approximating 
the  distribution  of  biweight-"t"  to  a  Student's  t  with  some 
degrees  of  treedom,  for  practical  purposes.  Hence,  we  shall 
make  the  correspondence 

(critical  point,  o)  --->  degrees  of  freedom  . 


When  we  examine  the  performance  of  biweight-" t"  on  samples 
of  only  five  observations  (Exhibit  1),  we  are  initially  dis¬ 
appointed  with  the  results.  Not  only  do  we  see  low  efficiencies 
in  the  lengths  of  the  confidence  intervals,  but  the  matched 
degrees  of  treedom  are  unusually  low.  It  appears  that  the 
numerator  ha*  extremely  heavy  tails;  hence,  is  matched 

to  a  Student's  t  with  few  degrees  of  freedom. 


3 ■  "Unusual"  Gaussian  samples . 

In  Exhibit  2  we  consider  the  (swindled)  estimate  of  one 
tail  probability.  Notice  that  these  samples  have  been  sin¬ 


gled  out  because  they  resulted  in  unusually  large  estimates 
All  of  these  samples  have  the  property  that  three  of  the 
five  observations  (just  over  half)  are  extremely  close 
together,  with  the  other  two  being  far  enough  away  that  the 
bisguare  function  assigns  them  zero  weight.  Such  Gaussian 
samples,  although  moderately  rare,  do  occur  with  mace  than 
21  frequency.  In  these  cases,  is  sure  to  grossly 

underestimate  <x ,  and  <<  Var  (numerator ) ,  since  the 
blsguare  operates  as  for  n«3  with  extremely  small  variance. 
Any  reasonable  robust  estimate  of  scale  would  perform  like¬ 
wise.  Exhibit  2(b)  presents  location  and  scale  estimates 

for  more  "typical”  Gaussian  samples.  In  these  samples, 

2 

is  much  closer  In  value  to  the  usual  sample  s  j  hence,  good 
performance  in  biweight-*t"  is  expected. 

4.  Quant  tfying  the  behavior  of  "unusual*  samples. 

If  we  can  improve  the  estimate  of 
var (numerator ) 

in  these  problematic  Gaussian  samples,  we  may  hope  that  a 
similar  improvement  may  be  used  when  the  underlying  distri¬ 
bution  is  not  Gaussian.  We  therefore  need  a  measure  by 
which  to  classify  the  "unusual"  samples.  Returning  to  the 
formula  (9)  for  a^,  two  possibilities  for  such  a  aeasuce 

are  suggested  byt 
n 

a)  S  w(u4) 
t-1 


♦ 


-  9  - 


-  10  - 


Of 

n 

b>  3  +'(u,) 
i*l 

where 


We  are  particularly  interested  in  distinguishing  those 
samples  for  which  one  or  sore  of  the  observations  are  far 
from  the  estimated  center.  This  corresponds  to  |u^|  >  1, 
for  which  V  (u^) *w(u.) -0.  Oue  to  the  nonotonicity  of 
w  ( )  u I ) /  smaller  values  of  the  weight  function  always  indi¬ 
cate  Increasingly  greater  distance.  It  appears  likely  that 
n 

5  w  ( u  j )  will  be  a  note  informative  ancillary  statistic  than 
i*l 
n 

S  V  . 

i-L 

Exhibit  3  shows  stem-and-leaf  plots  for  the  values  of 
n 

5  w(u.)  for  the  three  sampling  situations.  The  unusual 
i-1 

Gaussian  saeplea  described  above  all  fall  aeong  the  saeplea 
n 

for  which  5  w(u.):2.94.  The  eajority  of  the  saeplea  have 
i-1 
n 

3  wiu .)  14.80,  for  which  3^  per  formed  adequately.  The  case 
1*) 

n 

where  3  v(Uj):j.8B  corresponds  to  one  observation  being 
t-l 

treated  essentially  as  an  outlier.  The  stes-and-leaf  plots 


suggest  that  the  three  cases  can  be  specified  in  teres  of  a 
range  of  one  unit  in  the  value  of 
n 

W  >  2  w(u.)  <1\ 

i-l 

Since  we  would  like  to  choose  the  Interval  so  as  to  most 
clearly  differentiate  aeong  these  samples,  we  choose  end¬ 
points  where  the  density  of  W  is  lows 

two  "false  outliers"*  W  <  3.3 

one  "false  outlier":  3.3  <  W  <  4.3  (8) 

no  "false  outliers":  W  >  4.3. 

Here,  "false"  alludes  to  the  fact  that  these  observations, 
although  some  distance  from  the  bulk  of  the  sample,  ate 
nonetheless  bonafide  observations  from  the  same  distribution 
as  the  others.  Por  the  case  of  One-Wild  where 
3.3  <  W  <.  4.3»  the  outlier  does  in  fact  usually  correspond 
to  the  wild  shot  ifrom  s  N(0,  100)  distribution).  Hen¬ 
ceforth,  it  will  be  convenient  to  analyse  our  results  for 
n-5  not  only  by  situation  but  by  slice.  A  si  ice  la  defined 
by: 

a)  n,  a  given  number  of  observations: 

b)  P,  a  dlstr ibutlonal  situation:  (9) 

c)  a  range  of  valuos,  w^  and  Wy,  for  which 

wL  <  W  <  wu. 

Por  a  more  detailed  analysis  of  the  effect  of  W  on  the 
biwe lght-"t"  distribution,  we  generated  nine  slices  of  600 
samples  each,  where 


« 


9 


-  11  - 


-  12  - 


a)  n-5  Finally,  the  large  discrepancy  among  situations  in 

b)  F  *  Gaussian,  One-Wild,  or  Slash  ave(S^)  for  the  low-weight  slices  suggests  that  a  deeper 

c)  W  <  1.3,  3.3  <  W  <  4.3,  W  >  4.3,  look  at  the  behavior  of  these  samples  is  required. 


Exhibit  4  tabulates  the  estimated  frequencies  for  each 

2 

slice,  and  the  average  values  of  the  biweight  and  based 
on  the  600  simulated  samples.  We  see  that  a  low-weight 
slice  for  n^-S  is  relatively  Infrequent,  occurring  in  21-5% 
of  all  samples  from  our  situations,  yet  the  frequency  is 
just  large  enough  to  produce  the  low  efficiencies  in  the 
biweight-"t"  intervals  of  Exhibit  l.  Panel  8  of  Exhibit  4 
reveals  that  indeed  the  use  of  the  biweight  in  the  numerator 
of  "t*^,  despite  its  deflated  scaling,  is  not  the  real 
problem,  as  its  variance,  even  in  the  low-weight  Gaussian 
samples,  is  only  slightly  more  than  twice  the  variance  of 
the  optimal  mean.  The  biweight  is  a  big  success  in  the 
high-weight  samples;  notice  that  the  variance  of  the  optimal 
mean  in  the  Gaussian  situation  is  nearly  attained,  and  that 
in  all  high-weight  slices, 

ave  (denominator  of  2-var  (numerator  of  "t"^)  (10) 

In  the  medium-weight  slice,  Bqn.  (10)  already  approximately 
holds  for  the  more  sttetched-tal led  distributions,  but  it  is 
off  by  nearly  a  factor  of  10  in  the  Gaussian  situation.  In 
order  to  achieve  correspondingly  good  results  for  all 
medium-weight  slices,  it  is  likely  that  we  will  need  to  be 
conservative  in  some  places. 


».  Dlgt  ess  ion ;  Gt  anul ar  tty  of  the  we  ight  d  istr ibut ion . 

It  is  worth  commenting  on  the  granularity  of  the  dls- 
n 

ic  ibut  ion  of  W  «  5  w(u^)  for  the  three  situations.  This 

i-i 

tendency  is  partly  due  to  ouc  scale  estimate, 

*  ■  3bt  •  nl/2-sbi 


w(ut)  »  w((*i“Tbi)/(6&>) 

is  a  rather  extreme  case,  consider  the  following  estimate  of 

h 

fr-  (1/6 ) *  in  (Ixj^Xjl, lx2-x3l, lx4-x3l, U5-x3l) 
vhere  the  sample  x  is  assumed  ordered  (*^<.*2-*  * ’-*5^ '  Then 


i  ince  w  (u  . }  -0  for  all  l  except  1-3,  when  w(u3)-i.  Hence, 


n 

W  •  5  w(Uj)  -  1  -  constant, 

l-l 

regardless  of  any  further  characteristics  of  the  sample. 


« 


-  13  - 

That  this  is  a  rather  silly  estimate  Car  cr  can  be  seen  from 
the  following  two  fabc icated  samples: 


a) 

-1  .6, 

-0.8, 

-0.6, 

0.4, 

1.0 

•»>  &  » 

0.01 

b) 

-0.9, 

-0.8, 

-0.6, 

0.4, 

0.4 

■»>  <y  a 

0.  01 

Nonetheless,  the  example  does  serve  to  indicate  that  the 
continuity  of  the  density  function  of  W  is  highly  dependent 
upon  choice  of  scale.  It  is  quite  possible  that  there 
exists  a  choice  of  scale  foe  which  W  has  a  somewhat  smoother 
density  function.  Pot  reasonably  efficient  estimates  of 
scale,  however,  its  density  is  likely  to  have  modes 
separated  roughly  by  one  unit  (on  the  weight  scale).  The 
cutoff  points  we  have  selected  in  Equation  (9)  are  likely  to 
be  satisfactory  (l.e.,  to  come  at  very  low  densities)  for 
the  weighting  based  on  any  reasonable  scale  estimate. 


-  14  - 

PART  0:  COMPARTMENTALIZING:  SLICES. 

6.  A  scaled  biwe  iqht-*f  for  al  tees. 

Since  out  three  weight  classes  in  each  situation 

vaguely  represent  the  degree  to  which  S*4  falls  as  an  esti- 

Ute  of  the  wet  lance  ot  the  biwelght,  a  scaled  version  of 

•t*  ,,  conditional  on  a  given  weight  slice,  sight  have  a 
bi 

distribution  which  la  sore  stellar  to  a  Student's  t.  That 
is,  we  would  like  to  find  a  scale  factor,  K,  such  that 

"t" 

P(  |j-~  >  a  I  wLtwrwu  ,  n  )  111) 

T 

"  -  a  1  wl.<Wiwl!  *  "  1  -  1,1  V  *  ‘  1 

where  both  *  and  IT  Bay  depend  on  W  -  Slwelghts)  and  on  the 
sample  size  n. 

One  choice  of  it  is  suggested  by  the  values  In  Exhibit 
4(b).  If  we  want  to  Insist  that  l^n-1,  and,  in  addition, 
that  (10)  hold  approximately  in  all  situations,  we  would 
choose  our  scale  factors  as  follows: 


conservative  K 

Gaussian 

One-Wild  Slash 

(max  of  three) 

low  W 

IS. 49 

8.4S  2.21 

IS. 49 

medium  W 

4.16 

1.05  0.90 

4.16 

high  W 

0.94 

0.82  1.07 

1.07 

While  these 

scale  fectors  are  all  of  the 

same  order  in 

the 

medium-  end 

hlgh-we ight 

slices,  clearly 

we  may  be  much 

too 

conser  vat  We 

in  the  low- 

-weight  slice.  Furthermore,  it 

is 

4 
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not  altogether  certain  that  a  matching  of 

ave (denominator 2)  ~  var (nuraer ator ) 

will  transform  a  blweight-*t"  distribution  Into  one  from  the 
Student'e-t  family. 

An  adaptive  alternative  may  be  based  on  dealing  with 
the  actual  critical  points  from  biwe lght-"t".  Ideally, 

»  constant  foe  all  4^ 

or,  equivalent ly, 

Y  j  ( tO  -  log(*t"<4i)/t^<c<il)  •  constant  for  all  of  4 
A  least  squares  approach  would  minimize 

»  2 

5  lY.Urt  -  constant  \  (12) 

i-1 

whence 


constant  »  avefy^dO)  •  y(W. 

Tnen,  minimizing  (12)  would  be  equivalent  to 

%n  .’l* 

2 

where  a  [{/)  is  our  sample  variance  formula 

2  q  2 

•  -  S  (y 4 C to -y < to )  /(9-1) 

1-1 

However,  our  y  ^  ( (/)  are  not  independent,  and,  even  if  they 
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were,  the  usual  sample  s^  would  be  vetoed  on  the  basis  of 

2 

its  overall  nonrobustness.  As  we  mentioned  in  15),  s^  per¬ 
forms  well  with  a  sufficient  number  of  observations  [here  we 
will  use  g«)0,  where  9  is  the  number  of  values  of  y  4  <  tO  J • 
Moreover,  Portnoy's  results  (|1H|)  indicate  that  using  the 
redescending  '(-function  may  help  us  with  at  least  one  type 
of  dependence  among  the  observations.  Let  us  therefore 
choose  {/q  by 

</Qi  s2biU/Q)  is  a  minimum.  (13) 

Secondly,  we  select  a  more  conservative  value  for  the  con¬ 
stant  by  ^-signing 

log (scale  factor)  •  |<t<gyi(lV  1141 

7.  Log (scale  factors) ,  b^  si  ice . 

Exhibit  5  summarizes  the  degrees  of  freedom  and  log 
(scale  factors)  for  the  slices,  both  for  and  for  Vq 

chosen  via  (13).  The  closeness  of  our  fit  to  a  Student 4s  t 
on  (S  degrees  of  freedom  may  be  viewed  graphically  by  plot- 
t  ing 

log  ((scaled  "t"  (c^)  )/t^Wj)  >  vs  -log(4[) 

with  one  standard  deviation  "confidence  limits"  obtained 
from  the  curves 
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loq((scaled  eccoc  J /t  ^(ctp  »  vs  -logical 

Notice  in  these  plots  (Exhibit  6)  that: 

O)  a  negative  value  of  the  ordinate  indicates  a  con* 
servative  fit,  and 

(2)  a  negative  slope  suggests  a  larger  value  of  (/. 

the  most  successful  approximation  is  the  high-welght 
slice,  for  which 

(b iwe ight -"t "I/O. 95  approx imatejy  distributed  as  t4 
(C.f.  Exhibit  6(a)). 

For  the  medium-weight  slice,  a  uniform  rescaling  of  the 
biwe  ight-'t"  distribution  aore  closely  Batches  a  Student's  t 
on  3  d.f.  One  Bight  argue  that  the  Gaussian  sample  shows  a 
"suspect*  outlier,  and  a  conservative,  albeit  wasteful, 
approach  is  to  allow  ourselves  one  fewer  degree  of  fceedoa. 
There  is  only  a  snail  probability  that  we  will  waste  this 
valid  observation  in  the  Gaussian  case  (0.020  from  Exhibit 
4).  Of  course,  there  is  a  much  greater  likelihood  of 
obtaining  a  med ium-we ight  One-Wild  sample?  in  this  case, 
inference  based  on  four  of  the  observations  is  a  sensible 
procedure  in  the  absence  of  knowledge  of  the  kind  of  contam- 
(nation.  Exhibit  ?(a)  shows  the  relative  improvement  in 
comparing  our  scaled  *t*  points  to  Student's  t  on  3  d.f. 

The  scale  factors  for  the  low-weight  slice,  however, 
are  still  radically  different.  Not  surprisingly,  "t"  needs 
to  be  adjusted  more  drastically  when  the  underlying  dlstrl- 
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bution  is  truly  Gaussian.  Tor  a  decent  matching  at  the  less 
extreme  tail  areas,  a  scaled  "t"  compared  to  2  d.f.  otters  a 
possible  approach  (Exhibit  7(b)),  but  the  approximation  is  still 
far  from  good. 

Comparing  the  graphs  for  the  two  low-weight  fits,  one  pos¬ 
sible  procedure  is 

for  .05  <  a  <  .0005  ,  compare  "t"/21.0  to  Student's  t2, 

for  .0005  <  a  <  .00001,  compare  "t"/91.2  to  Student's  l^. 

It  is,  however,  worthwhile  to  characterize  the  differences  in 
the  three  low-weight  classes.  One  difference  is  apparent  from 
Exhibits  5  through  7:  the  scale  factors  and  approximations  for 
One-Wild  and  Slash  are  very  similar  to  one  another  and  each  is  con¬ 
siderably  different  from  the  Gaussian.  If  we  had  a  method  whereby 
we  could  discriminate  the  Gaussian  samples  from  those  whose  under¬ 
lying  distribution  has  more  stretched  tails,  we  might  be  more 
successful  in  adaptively  scaling  biwe ight-nt" .  This  idea  is 
pursued  in  greater  detail  in  (10 3. 


rowed .  ) 


FART  C:  CLOSE. 

a.  Conclusions  for  n=5. 

The  initial  aim  of  this  study,  that  of  constructing  a  valid 
confidence  interval  for  the  center  of  a  population  for  which  we 
have  only  five  observations ,  led  to  a  more  ambitious  goal  of 
characterizing  the  distribution  of  the  biweight-"t"  statistic. 

The  case  with  r.'S  is  perhaps  the  most  difficult  of  all:  there 
are  too  many  observations  to  develop  an  analytic  solution,  yet 
so  few  that  the  likelihood  of  obtaining  potentially  misleading 
samples  cannot  be  ignored.  In  searching  for  a  more  complete  des¬ 
cription  of  the  tail  behavior  of  "tM  on  five  observations,  we 
discovered  that  a  characterization  based  on  the  sum  of  the  (bi¬ 
weight)  weights  offers  a  more  satisfying  approach. 

When  we  can  borrow  information  on  width  from  several  samples 
2 

we  may  compute  both  T^  and  using  a  scale  estimate  pooled 

from  all  samples: 


s  ,  »  (9s  »\|Jn’ 
pool  p 

jn 

5  +  (u 
t-1 

i> 

Jn 

Jn 

f  5  + 

i »] 

5  +’  (Uj)) 
l-l 

Sp  »  pooled  MAD 

u  * 

In  this,  case,  as  Exhibit  8  reveals,  biweight-"tl*  performs  fairly 
well.  (This  Exhibit  tabulates  only  the  matched  degrees  of  free¬ 
dom  and  ECU.  efficiencies  for  the  sake  of  brevity;  J  refers  to 
the  number  of  samples  from  which  scale  information  has  been  bor- 


It  is  not  always  clear,  however,  if  and  when  we  may  bor¬ 
row.  An  unjustified  usage  of  borrowing  may  be  extremely  mislead¬ 
ing.  In  this  regal'd,  the  sum  of  the  weights  may  lend  insight: 
an  unexpectedly  low  value  of  W  may  caution  us  to  treat  this  sample 
separately  from  the  rest  and  not  use  it  in  borrowing  width  infor¬ 
mation.  The  borrowing  issue  plays  a  more  important  role  in  the 
two-sample  problem  [10]. 

In  the  absence  of  additional  samples  of,  say,  five  observa¬ 
tions,  or  any  additional  width  information  for  our  sample  at  hand, 
a  conservative  approach  to  the  interval  estimation  problem  for  the 
single  sample  of  five  observations  would  be 
5 

A)  For  Z  w(u^)  >  4.3,  use  t 4 (a) * (0 . 95‘S^. ) 

for  the  allowance  in  the  confidence  interval. 

As  we  saw  in  Section  6,  the  resulting  confidence 

interval  performs  very  well; 

B)  For  3.3  <  Ew(u^)  <  4.3,  use  either  (a ) * < 7 . 1 ’Sfa . ) 

or  t , (a) * (4 . 8 *S. . )  for  the  allowance; 

3  bi 

C)  For  Ew(Uj )  <  3.3, 

(i)  for  .05  <  a  <  .001  ,  compare  Ht"/21.0  to  Student's  t? 
for  .001  <  a  £  . 00001,  compare  "t"/91.2  to  Student's  t^ 

(ii)  consider  the  improvements,  based  on  additional  an¬ 
cillary  statistics,  given  in  [10); 

(iii)  pray  for  more  information. 


Tall  Pr. 

Crit.  Pt. 

Stnd.  Error 

D.F. 

ECIL 

Efficiency 

0.00001 

222.6 

{  8.206) 

2.0 

187.3 

1.10 

0.000025 

219.6 

(  7.654) 

2.0 

184.8 

0.71 

0.00005 

206.7 

( .7.596) 

2.0 

173.9 

0.56 

Dist'n. 

0.00010 

186.0  | 

[  7.397] 

1.9 

156.5 

0.47 

Gaussian 

0.00050 

118.1  1 

5.897] 

1.7 

99.34 

0.53 

0.00100 

86.79  (  4.087] 

1.6 

73.03 

0.68 

0.00500 

26.93  (  1.098] 

1.4 

22.66 

2.92 

0.01000 

13.20  (  0.479] 

1.6 

11.11 

8.04 

0.02500 

4.325  (  0.126] 

2.0 

3.639 

41.16 

0.05000 

2.650  (  0.061] 

2.3 

2.230 

64.64 

0.00001 

133.8  1 

[  4.843)  2.6 

169.9 

1.33 

0.000025 

120.5  l 

4.804)  2.2 

153.0 

1.04 

0.00005 

109.0  I 

;  4.702)  2.0 

138.4 

0.89 

Dist'n. 

0.00010 

95.77  1 

;  3.853)  2.0 

121.6 

0.81 

One-Wild 

0.00050 

62.60  1 

[  2.633 

1  1.9 

79.47 

0.83 

0.00100 

47.70  1 

2.231] 

1  1.8 

60.55 

0.99 

0.00500 

16.10  1 

[  0.626)  1.7 

20.44 

3.59 

0.01000 

8.170  < 

;  0.281 

1  1.8 

10.37 

9.22 

0.02500 

3.390  1 

0.092, 

»  2.7 

4.304 

29.43 

0.05000 

2.196  1 

;  0.047)  3.6 

2.788 

41.36 

0.00001 

294.188 

[13.037] 

2.0 

763.308 

5.70 

0.000025 

258.162 

[12.597] 

2.0 

669.833 

5.78 

0.00005 

227.870 

[12.035 

1.9 

*  591.238 

7.18 

Dist'n. 

0.00010 

210.231 

[10.899] 

1 .9 

545.469 

6.66 

Slash 

0.00050 

102.154 

,  6.298] 

1.7 

265.051 

5.44 

0.00100 

53.148 

[  3.305] 

1.8 

137.901 

5.88 

0.00500 

15.011 

[  0.553 

1.7 

38.947 

7.42 

0.01000 

8.820  < 

[  0.281] 

1.8 

22.886 

9.17 

0.02500 

3.954 

[  0.106] 

2.2 

10.258 

11.58 

0.05000 

2.472 

[  0.056] 

2.7 

6.415 

10.06 

Exhibit  1:  Results  on  one-sample  biweight  -  "t"  , 


n  *  5 


Exhibit  2 


(A)  p{"t"b.  >  3.747}  for  14  "unusual"  Gaussian  samples. 


Sample 

Number 

x(l) 

x(2) 

x(3) 

x(4) 

x(5) 

Tbi 

Sbi 

X 

^sample 

s2  , 

sample 

P 

25 

-0.552 

-0.434 

-0.479 

0.227 

0.658 

-0.504 

0.024 

-0.126 

0.242 

100.628 

0.44  37 

59 

-1.338 

-1.292 

-0.012 

0.010 

0.123 

0.040 

0.043 

-0.502 

0.333 

59.957 

0.4303 

159 

-0.093 

-0.088 

-0.070 

0.132 

1.667 

-0.084 

0.007 

0.310 

0.342 

2439.455 

0.484  3 

165 

-0.958 

-0.940 

-0.930 

-0.188 

2.111 

-0.942 

0.008 

-0.181 

0.591 

5041  .095 

0.4935 

178 

-1.800 

-1.349 

-0.056 

-0.030 

0.008 

-0.026 

0.019 

-0.645 

0.386 

409.504 

0.4755 

299 

-1.902 

-0.540 

0.540 

0.551 

0.607 

0.566 

0.021 

-0.149 

0.489 

521.027 

0.4799 

328 

-0.994 

-0.816 

-0.773 

1.312 

2.012 

-0.860 

0.070 

0.143 

0.629 

81.376 

0.4775 

388 

-0.480 

-0.460 

-0.448 

0.608 

1.400 

-0.463 

0.010 

0.124 

0.381 

1527.198 

0.4362 

444 

-0.270 

-0.173 

-0.158 

0.834 

1.014 

-0.200 

0.036 

0.249 

0.277 

58.971 

0.4271 

511 

-0.513 

0.526 

0.554 

0.556 

1.026 

0.546 

0.010 

0.430 

0.254 

645.000 

0.4678 

515 

-1.454 

-1.442 

-1.315 

-0.115 

0.765 

-1.404 

0.045 

-0.712 

0.446 

98.448 

0.4531 

535 

-0.898 

-0.879 

-0.775 

0.396 

1.187 

-0.851 

0.039 

-0.194 

0.422 

117.617 

0.4549 

575 

-0.265 

0.517 

0.544 

0.552 

2.094 

0.533 

0.011 

0.688 

0.384 

1210.235 

0.4652 

604 

-2.877 

-2.346 

-0.091 

-0.013 

0.039 

-0.022 

0.038 

-1.058 

0.640 

278.543 

0.4906 

mean 

-0.262 

0.027 

-0.116 

..  0.415 

899.220 

0.4669 

std  err (mean) 

0.162 

0.005 

0.130 

0.035 

360.291 

0.0060 

Exhibit  2  (continued) 


(B)  pCt" 

'  >3.747} 

bl 

for  12 

"typical" 

Gaussian  samples. 

sample 

number 

x(l) 

x(2) 

x(3) 

x(4) 

x(5) 

Tb1 

Sb1 

x 

Ssample 

s2  . 

sample 

Sb1 

P 

606 

-0.657 

-0.022 

0.217 

0.345 

2.733 

0.131 

0.404 

0.523 

0.579 

2.052 

0.0030 

616 

-0.472 

0.361 

0.615 

0.693 

0.741 

0.568 

0.137 

0.388 

0.225 

2.679 

0.1440 

6T7 

-0.321 

0.234 

0.247 

0.703 

0.892 

0.356 

0.224 

0.351 

0.211 

0.886 

0.0301 

618 

-0.457 

0.633 

0.852 

1.070 

1.114 

0.700 

0.293 

0.640 

0.290 

0.979 

0.0075 

619 

-0.633 

-0.178 

0.247 

0.395 

0.871 

0.143 

0.272 

0.141 

0.256 

0.887 

0.0114 

621 

-0.178 

0.183 

0.366 

0.575 

1.034 

0.393 

0.211 

0.396 

0.202 

0.912 

0.0333 

622 

-1.259 

-0.574 

0.129 

0.612 

0.912 

-0.025 

0.424 

-0.036 

0.395 

0.870 

0.0002 

623 

-1.273 

-0.726 

-0.686 

-0.285 

0.151 

-0.565 

0.251 

-0.564 

0.238 

0.897 

0.0177 

624 

-0.771 

-0.507 

-0.406 

0.095 

0.983 

-0.143 

0.322 

-0.121 

0.310 

0.923 

0.0035 

627 

-0.086 

0.000 

0.289 

0.455 

1.091 

0.336 

0.219 

0.350 

0.210 

0.917 

0.0333 

629 

-0.539 

-0.538 

-0.045 

0.715 

1.169 

0.143 

0.367 

0.152 

0.342 

0.868 

0.0010 

631 

-1.417 

-0.996 

-0.378 

0.130 

0.192 

-0.487 

0.340 

-0.494 

0.314 

0.855 

0.0022 

mean 

0.129 

0.289 

0.144 

0.298 

1.144 

0.0240 

std  err 

(mean) 

0.111 

0.025 

0.110 

0.031 

0.169 

0.0120 
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Sxhibit  3 

Stem-and-leaf  plots  foe  W=5w(u.) 
foe  thcee  sampling  situations  (n-5) 


(i)  640  Gaussian  samples: 


29 

30 

31 

32 

33 

34 

35 

36 

37 
33 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 


4444444444444444444 


3888883388333 

0677 

15 

347 

07 

136 

56 

57 

01456777773338888339999999999999999999999999999999+15 

00000000000000000000000000000000000000000000003000+99+ 

01244 


(ii)  640  One-Wild  samples: 

291444444444444444446 

30  I 

31  ! 

32! 

331 

34124 

351 

36  1 

37  1 

38  I  36633838338338838888838888  388833  333338338883898888+99+ 
3910122344555799 

40 | 3446 
41  I  3579 
421 

4310246 

441389 

4512235 

4611133357338 

47  1011  2234  5  556777777738883399999999999999999999999999+35 
43  I  00000000000000000000000000000000000000000000000000  +  99  + 
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Exhibit  3  (cont.) 
Stem-and-leaf  plots  foe  W=2w(u.) 
foe  thcee  sampling  situations  (n-5) 


(iii)  640  Slash  samples: 


291444444444444444444444444444444444446 

30  I  9 

31  I 
321 
331 
341 
351 
36  I 
3713 


38  I  5778393888889983898338888  83899938338888  389939983333+91 

391002235539 

4011233473 

41  I  56 

42  I  5 

43  1 

44  I  3459 
45135579 
4613839 

47 111233455566777839339999999999999999999999999999999+22 
48 ) 00000000000000000000000000000000000000000000000000+99+ 
49  I  4 


( iv)  Numbec  of  samples  in  thcee  weight  gcoups: 


oist'n 

1  W< 3 . 3 

1  3.3  <W< 4 . 3  1 

W>  4 . 

Gaussian 

1  19 

1  23  1 

599 

One-Wild 

1  13 

1  323  I 

299 

Slash 

1  37 

1  151  I 

442 
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Exhibit  4 

Information  of  slices  foe  n  =  5 


M  Relative  frequencies  (in  %)  of  slices 
(standard  errors  in  parentheses) 


Gauss  ian 

One-Wild 

Slash 

W<  3 . 3 

2.43 

(0.10) 

2.30 
(0.11  ) 

4.57 

(0.13) 

3. 3<W<4. 3 

2.76 

(0.10) 

45.57 

(0.34) 

21  .1  4 
(0.36) 

W>  4 . 3 

94.81 

(0.14) 

51  .63 
(0.34) 

74.29 

(0.33) 

3)  Some  summacy  values  on  Tfai  and  S^,  by  slice 


W<3 . 3 


3 . 3  <W< 4 . 3 


Gauss ian 

One-Wild 

Slash 

vac  <Tm) 

0.480 

(.014) 

0.500 
(.011  ) 

1  .272 
(.065) 

0.002 
(.0001  ) 

0.007 
(.0001  ) 

0.260 

(.032) 

vac<Tbl) 

0.380 
(.  005) 

0.253 
(.001  ) 

2.242 

(.220) 

*«<s  bi1 

0.020 

(.015) 

0.230 

(.002) 

2.77] 

(.306) 

vac(Tbil 

0.203 

(.004) 

0.323 

(.033) 

9.450 

(5.050) 

ava(S^) 

0.223 

(.003) 

3  .219 
(.052) 

9.300 

(2.647) 

W>  4 . 3 
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Exhibit  5 

Log (scale  factors),  by  slice 


h)  Fitted  log(scale  factors)  (=log(K))  and 
degrees  of  freedom  (  =  (/): 


Gauss  ian 

One-Wild 

Slash 

Low 

</ 

2 

2 

3 

We  ight 

log  (K) 

1.59 

1.22 

j 

1.53 

Med ium 

</ 

3 

4 

3 

We  ight 

log (K) 

0.68 

0.23 

0.17 

H  igh 

</ 

5 

5 

4 

We igh  t 

log (K) 

0.064 

-0.035 

-0.10 

3)  Log 

(scale  factors)  for  degrees  of  freedom  =  4 

Gauss ian 

One-Wild 

Slash 

Low 

(/=4: 

1.79 

1.96 

1.71 

Med ium 

(/=4 : 

0.35 

0.28 

0.34 

H  igh 

(/=4: 

-0.016 

-0.12 

-0.10 

C)  Log 

(scale  factors)  for  (/* 

-]+[-! +5we ights] : 

Gauss ian 

One-Wild 

Slash 

Low 

V-2: 

1  .  32 

1  .22 

1.19 

Med ium 

(/=3: 

0.69 

0.13 

0.13 

H  igh 

(/=4: 

-0.016 

-0.12 

-0.10 

on  medium-weight  slices,  n=5. 
(g=Gaussian;  w=0ne-Wild;  s=Slash; 
including  1  std.  dev.  in  "t" (a) ) 


t"(a)/91. 


t" (a)/4  . 
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on  medium-weight  slices,  n=5. 

(g  =  Gaussian w=0ne-Wild;  s  =  Slash; 
including  1  std.  dev.  in  "t" (a)) 


t" (u)/21. 
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Exh  ibit  3 

Giweight-"t"  on  n=5:  3or cowed  numerators 
and  denominators,  using  pooled  ?i3YMV. 


»boc rowed  | _ 2 _ | _ 3 _ | _ 4 _ I _ 5 


Id 

.  f . 

eff . 

I  d.f. 

eff. 

Id.  f . 

eff . 

Id.f . 

ef  f . 

Gauss  ian 

.0000] 

1 

5.2 

97.0 

|  3.5 

32.0 

|  9.0 

35.1 

1 1  9 . 1 

36.3 

.00005 

1 

5.2 

96.2 

|  3.4 

34.3 

1  9  .  5  ' 

88.6 

1 1  3.5 

35.7 

.000] 

1 

5.0 

95.1 

|  3.5 

36.5 

1  9.3 

91  .1 

1 1  3.5 

35.3 

.0005 

1 

5.9 

98.4 

|  9.5 

90.8 

1 1  1  .7 

92.0 

1 1  9.3 

87.7 

.001 

1 

5.5 

99.3 

|  10.2 

99.0 

1 1  3 . 1 

93.5 

|  20. 3 

89.1 

.005 

1 

5.6 

]  0]  .6 

|  11.6 

101.1 

1 1  6 . 2 

107.6 

|  26.3 

93.2 

.0] 

1 

5.9 

115.7 

|  12.5 

103.7 

|1  3.4 

110.4 

|  33.7 

95.0 

.025 

1 

6.5 

127.4 

|  14.5 

103.7 

|25.1 

113.7 

|  32- 0 

97.5 

.05 

1 

7.0 

1  33.4 

1  13.4 

105.7 

147.4 

115.2 

i  CO 

99.4 

One-Wild 


00001 

|  5.2 

58.1 

1 

8.9 

70.9 

1  9.3 

92.0 

|  24 . 4 

35.1 

00005 

1  5.3 

60.5 

1 

9.0 

70.1 

1  9.3 

91  .5 

|23.0 

34.4 

0001 

|  6.0 

63.1 

1 

9.0 

56.5 

1  9.7 

91 .1 

|30.4 

34.2 

0005 

|  5.2 

54.0 

1 

3.3 

64.0 

|13.4 

90.2 

|  4  2 . 3 

33.3 

001 

1  5.3 

68.1 

1 

3.3 

62.4 

|16.6 

94.7 

|55.3 

83.3 

005 

|  5.3 

75.2 

1 

12.9 

69.2 

|34.1 

96.1 

1  o° 

33.2 

01 

|  6.6 

31  .2 

1 

16.7 

71  .0 

(76.0 

96.9 

1  oo 

32.8 

025 

|  3.5 

33.5 

1 

32.7 

71  .5 

00 

95.6 

1  oo 

32.1 

05 

111.5 

33.1 

1 

00 

70.9 

1  CO 

95.2 

1  00 

30.3 

Slash 


0000] 

|  7.3 

41.2 

(  10.0 

43.3 

1 1  2 . 2 

50.1 

|11  .2 

49.3 

0000  5 

1  7*7 

40.4 

|  10.0 

45.4 

[12.1 

52.2 

1 3  0 .9 

43.3 

900] 

1  7  •  5 

40.0 

|  9.3 

46.  .3 

|11  .3 

51  .1 

|11  .0 

43.2 

0005 

1  7  •  5 

42.1 

|  10.4 

49.9 

|11. 5 

49.2 

(11.7 

47.5 

001 

|  3.0 

44.2 

|  10.6 

52.2 

1 1  2.0 

43.0 

|  1  3  .9 

48.5 

005 

|  3.2 

50.0 

|  1 1  .1 

54.3 

1 1  2 . 4 

49.9 

|  12.1 

49.1 

01 

|  3.4 

51  .5 

|  1  1  .2 

53.1 

(15.4 

51  .3 

(17.5 

50.3 

025 

|U  .5 

52.3 

|  19.0 

60.7 

|  54.1 

52.2 

1  00 

51  .9 

05 

123.1 

52.3 

1  202.4 

63.6 

I  00 

55.5 

1  00 

53.9 
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